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ABSTRACT 



The fine-structure spacing de,{n) = z/^,„ — z/^+2,n-i for low-degree solar p 
modes of angular degree i and radial order n is sensitive to conditions in the 
deep radiative interior of the Sun. Here we present fine-structure spacings 
derived from the analysis of nearly five years of helioseismological data collected 
between 1991 July and 1996 February by the Birmingham Solar-Oscillations 
Network (BiSON). These data cover 9 < n < 28 for do{n), and 11 < n < 27 
for di{n). The measured spacings are much more precise and cover a greater 
range, than earlier measurements from BiSON data (Elsworth et al. 1990a). 
The predicted fine-structure spacings for a "standard" solar model are clearly 
excluded by the BiSON data (at ^ lOa); models that include helium and heavy 
element settling provide a much better match to the observed spacings (see also 
Elsworth et al. 1995). Since the inclusion of core settling in solar models will 
tend to slightly increase the predicted neutrino flux, the BiSON fine-structure 
data appear to reinforce previous conclusions, i.e., an astrophysical solution to 
the solar neutrino problem seems unlikely. 

Subject headings: Sun: interior - Sun: oscillations - Stars: evolution 
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1. INTRODUCTION 

The solar p-mode oscillations are manifestations of trapped standing acoustic waves 
in the solar interior. The low-degree (low-£) modes, which are globally coherent, possess 
radial wave functions that penetrate into the core of the Sun. Fine-structure spacings 
de{n) = vi^n — vi+2,n~i (where n is the modal radial order), constructed from the low-degree 
modes, are extremely sensitive to conditions in the deep solar interior. The use of these 
spacings to constrain models and theory is particularly appealing, since one avoids relying 
upon inadequate modelling of the outer solar layers. 

Here we present fine-structure spacings (io(n) and di{n) derived from the analysis 
of several years of solar Doppler velocity data collected between 1991 and 1996 by the 
Birmingham Solar-Oscillations Network (BiSON). Elsworth et al. (1990a), with pre-1990 
vintage BiSON data, demonstrated that the "standard" solar model provided a good 
match to the observed spacings. Further, a variety of WIMP models (e.g.. Cox, Guzik & 
Raby 1990) were shown to be in significant conflict with the observations. Elsworth et al. 
(1990a) therefore concluded that the helioseismological data provided evidence for the solar 
neutrino problem requiring a solution in particle - and not astro - physics, e.g., via mixing 
between different species of neutrinos or spin precession. The substantially superior quality 
of the BiSON data presented in this paper has allowed us to make even more precise - and 
demanding - comparisons with solar models. We demonstrate that models including helium 
and heavy element settling are in much better agreement with the new BiSON data (see 
also Elsworth et al. 1995). 
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2. DATA AND ANALYSIS 

We have used seven 8 month-long, and one 32 month-long time series, all constructed 
from BiSON data collected between 1991 July and 1996 February, for the analyses discussed 
in this paper. The modes have been fitted in the frequency domain as rotationally split, 
Lorentzian multiplets with an associated fiat background offset. The relative amplitudes 
of different m components within individual i multiplets were fixed with the whole-disc 
sensitivities of Christensen-Dalsgaard (1989). We have constrained the heights of peak 
components with the same \m\ to be equal within each fitted multiplet. A Levenberg- 
Marquardt technique was used to perform the fitting, minimizing a maximum likelihood 
function that refiects the x^ 2-d.o.f. statistics of the frequency spectrum (Chaplin et al. 
1996). Dziembowski & Goode (1996) indicate that m-dependent asymmetries in the real 
mode-multiplet structures might introduce systematic errors. Our mode-fitting procedure is 
most heavily weighted toward the sectoral components, which appear as the strongest peaks 
in our data. The calculations of these authors suggest that - considering, for example, 
doln) - the frequency differences between the i = and i = 2, m = ±2 components should 
not be strongly affected over the solar cycle. 

Fine structure spacings do{n) and di{n) were determined for each spectrum from the 
fitted frequencies. Formal errors on the spacings were derived from the formal uncertainties 
on the mode frequencies returned by the fitting program. Mean sets of spacings were 
constructed from the seven 8-month spectra by computing weighted averages according to 

N 



Y,di{n)i/a'l 



di{n) = ^^ , (1) 



El/ 



a? 



i=l 

where A^ is the number of independent measures at each n, i.e., N < 7, and o"j are the errors 
associated with each spacing di{n)i. Estimates of the external error on each weighted mean 
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[e.g., see Topping 1962) were calculated according to 



4xtK(n)]=t[iV-l]- 



\ 



^ i , (2) 



{N-1)J2 l/< 



i=l 



where, because of the small number of samples at each {£,n), a suitable correction factor, 
t[N — 1], appropriate to the desired la significance interval, has been drawn from the t 
distribution in order to correctly assess the confidence limits bracing the mean (e.g., see 
Book 1978, p. 85). Constructing the errors in the manner shown above reflects the scatter 
in the observables, i.e., an analysis indicates that the external errors are a factor of ~ 1.5 
larger than the internal errors, while preserving the additional information provided by 
the formal fitting uncertainties. A small number of data were omitted from a few of the 
mean spacing calculations on the basis of several "outlier" rejection criteria. We employed 
variants of the standard discordancy test (internally studentized extreme deviation from the 
sample mean; e.g., see Barnett & Lewis 1984, p. 167), omitting each measure in turn from 
the sample in the same spirit as the "jackknife" technique and also taking into account the 
internal weights in one of the tests. In addition, we also used the Dixon Q family of tests 
(Dixon 1951) in order to check for the presence of single (Dixon's rio) and double (Dixon's 
r2o and r2i) outliers. 

Final sets of fine-structure spacings were then constructed by combining the 32-month, 
and mean 8-month data: for (io(^), 32-month spacings were used for 9 < n < 13, and 
mean 8- month spacings for lA < n < 28; for di{n), 32- month spacings were used for 
11 < n < 13, and mean 8-month spacings for 14 < n < 27. The errors used for n < 13 
are the formal uncertainties from the fitted 32-month spectrum, while those above are the 
external errors from the 8- month averages (Equation ^. The averaged data sets can be 



found at [http : //bison . ph . bham . ac . uk 
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3. RESULTS AND DISCUSSION 

We have followed the convention of Elsworth et al. (1990a), and parameterized each 
set of data in terms of a straight line according to 

de{n) = Co + ci ■ {n - 21) , (3) 

where cq and ci are the coefficients of the fit, suitably normalized to the radial order datum 
n = 21 on the abscissa. We have performed the fit over the range 15 < n < 27 in order 
to facilitate a proper comparison with Elsworth et al. (1990a). The best-fit coefficients 
for the Elsworth et al. (1990a) and new BiSON data are given in Table The new data 
are represented graphically in Fig. |I|: here, the residuals, generated by subtracting the 
appropriate best straight-line fit from the mean fine-structure spacings do{n) and di{n), 
have been plotted. (See later for explanation and discussion of the solid, dot-dashed and 
dotted-line model-generated predictions.) The uncertainties on the new best-fit values are, 
as expected, substantially superior to those in Elsworth et al. (1990a), owing to the higher 
quality and quantity of the more recent data. The fitted cq coefficients for di{n) differ by 
some 3a (combined error): this is because the £ = 3 frequencies - the mode-peaks of which 
are substantially weaker than those at £ = 0, 1 and 2 - were less- well determined in the old 
BiSON spectra. 

Fitting statistics indicate that a straight line is an inadequate representation of both 
sets of data, i.e., with reference to Fig. |l|, there is clearly rather more structure present 
in each plot. The use of higher-order terms, or other more-complicated functions, in any 
parameterized description of do{n) and di{n) is therefore now required, given the accuracy 
of the measured spacings. 

We have searched for solar-cycle effects in both (io(n) and di{n). No significant 
variations were found. For do{n), we find that a variation in the mean spacing, as weighted 
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by the global view, over the range 14 < n < 26, is excluded at ~ 0.08 /xHz {3a); and for 
di{n), over the same range in n, it is excluded at ~ 0.09 /xHz {3a). These values should 
be compared with the corresponding ~ 0.45 /iHz change observed in the eigenfrequencies 
themselves (Elsworth et al. 1990b, 1994; Libbrecht & Woodard 1991; Regulo et al. 1994; 
see also Evans & Roberts 1992). From an asymptotic description of the acoustic modes 
(e.g., Tassoul 1980) it can be shown that (for n » t) the fine spacings are proportional 
to the gradient of the sound speed. Since n/ 1 is quite small for some of the observed 
low-degree modes that make up the presented fine-structure spacings, care must be taken 
in interpreting these data through the use of asymptotics. This assumed, the solar-cycle 
constraints on dQ{n) and di{n) imposed by the BiSON data - in each case an upper limit to 
any change of less than 1 per cent - would appear to place similar tight constraints on any 
implied variation of the sound-speed gradient through the deep solar interior over the solar 
cycle. 

We have compared the observed BiSON fine-structure spacings with those derived 
from three, slightly different solar models. These are referred to as models A, B and C 
in Table |l[ Model A is a "standard" solar model (model 1 from Christensen-Dalsgaard, 
Profiitt & Thompson 1993; official designation 14b. 14), which neglects settling effects, and 
was constructed with the CEFF equation of state (Christensen-Dalsgaard & Dappen 1992) 
and the Livermore (OPAL) opacity tables (Iglesias, Rogers & Wilson 1992). Model B 
(model 4 from Christensen-Dalsgaard, Proffitt & Thompson; official designation 14b. d. 20) 
again employs the CEFF equation of state and the OPAL opacities, in addition to helium 
and heavy-element settling with turbulent mixing. Model C (Reference model 'S' from 
Christensen-Dalsgaard, Dappen et al. (1996), which appears as model OPAL 1 in Basu et 
al. (1996); official designation 15bi.d.l5) incorporates settling effects - neglecting turbulence 
- and employs the OPAL opacities, in addition to the OPAL equation of state (Rogers, 
Swensen & Iglesias 1996). 



The fine-structure data for each model were fitted, as per equation ^, over the range 
15 < n < 27. An inspection of the fitted Cq coefficients hsted in Table |l] clearly indicates 
that model A is excluded at a high level of significance. While models B and C provide a 
better match to the observed data, the do{n) comparison favours model B, while that for 
di{n) favours model C The fitted slopes for the models are all fairly consistent: they lie 
roughly 3 and 1.5a from the new BiSON values for dQ{n) and di{n) respectively. 

The use of the fitted BiSON and model coefficients in assessing the relative merits 
of the models is, to some extent, misleading, i.e., the comparison is ultimately dependent 
upon the suitability of the parameterized description - here, a simple straight line - used to 
describe the trend in the spacings. As previously noted, the plots in Fig. ^ clearly show that 
a straight-line representation of the measured spacings is inadequate. We have therefore 
also employed a direct one-to-one comparison between the BiSON and model spacings, 
as indicated in Table |I[ For each model under consideration, fine-spacing residuals were 
computed at each n, in the sense of BiSON minus model. Mean, weighted residuals were 
then computed for a variety of ranges in n - those for the range 9 < ra < 28 for do{n), 
and 11 < n < 27 for di{n), are shown in Table |l]. The uncertainties on each fine-structure 
datum were used to weight the computation of the mean, and its error determined from the 
weighted scatter of the residuals. 

With reference to Table |l], the difference residuals imply that the spacings for the 
"standard" model (A) are significantly larger than the BiSON values. The mean residuals 
for do{n) and di{n) exclude model A at the 13 and ~ 9a levels. The models incorporating 
settling provide a much better match to the observed spacings; however, marginally 
significant differences do remain. The comparative do{n) statistics imply that model B is in 
slightly better agreement with the BiSON spacings than model C Fig. |1| appears to show 
that model B (dot-dashed line) provides a better match to the the BiSON do{n) data for 
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n > 17, while at lower n models B and C are comparable. (The weighted difference residuals 
for 9 < n < 17 are significant at the ~ 3 and ~ Aa levels for models B and C respectively; 
while for 18 < n < 28 they are significant at the ~ 0.1 and ~ 5a levels respectively.) For 
di{n), the difference residuals for models B and C are similar; model C is perhaps more in 
accordance with the observed values. 

As an aside, we note that the inclusion of turbulence in model B has little effect on 
the spacings. The introduction of turbulence in the model changes the sound-speed profile, 
between approximately 0.5 < r/R^ < 0.7, relative to a model that neglects its effects. 
However, the sound-speed profile - and the hydrogen and helium abundances - are very 
similar in the deep radiative interior. Christensen-Dalsgaard, Proffitt & Thompson also 
constructed a settling model that neglected turbulent effects: this model produces cq, ci, 
and mean residual values in good agreement with model B. 

The plots in Fig. |l|, and the comparative residual statistics, show clearly that model 
A is excluded - at a very high level of significance - by the BiSON data. While neither 
of the "settling" models considered provides a satisfactory match over the whole range in 
n for which the observed spacings are available, they do nevertheless match the BiSON 
data far better. The inclusion in solar models of non standard settling effects in the core 
will tend to raise the opacity of the deep radiative interior and increase the solar neutrino 
flux, albeit rather modestly (Christensen-Dalsgaard 1996). The results in this paper would 
therefore appear to reinforce the conclusions of Elsworth et al. (1990a), in arguing against 
an astrophysical solution to the solar neutrino problem. 

Christensen-Dalsgaard (1996) does, however, point out that certain models can be 
constructed (e.g., Antia & Chitre 1995) that match observed solar neutrino capture rates 
and the helioseismological data, but at the cost of employing what appear to be unrealistic 
assumptions regarding the physics of the solar interior, for example, models that arbitrarily 
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"juggle" the comparative contributions of a reduced core opacity (which will tend to reduce 
the fine-structure spacings) and mixing in the core (which will tend to increase them). 
At present there appears to be little justification for constructing such models. We add 
an additional note of caution by indicating that the characteristics of the p modes in 
themselves - which depend upon the mechanical properties of the solar interior - cannot 
uniquely define the internal solar temperature and, by implication, the expected neutrino 
fiux. Nevertheless, helioseismological data, such as those presented in this paper, reinforce 
our previous conclusions, i.e., an astrophysical solution to the solar neutrino problem seems 
to be excluded. If the various nuclear cross sections are correct, then electron neutrinos 
would therefore seem to disappear between the solar core and terrestrial detectors. Possible 
causes are: oscillations between various species, or the precession of a magnetic moment, 
both of which require neutrinos to have mass (e.g., see Bahcall 1988). 

We would like to thank all those who are, or have been, associated with the BiSON 
global network. In particular: H. B. van der Raay, H. Williams, J. Litherland and R. Lines 
in Birmingham and P. Fourie at SAAO: and our hosts R. Stobie (SAAO); the Carnegie 
Institute of Washington; the Australia Telescope National Facility (CSIRO); E. J. Rhodes 
(Mt. Wilson, California); and members (past and present) of the lAC, Tenerife. We are 
indebted to J. Christensen-Dalsgaard for providing the model-calculated frequencies. Some 
of the research presented above made use of the Birmingham node of STARLINK. The 
10.7-cm radio fiux data were obtained from the World Data Centre. BiSON is funded by 
the UK Particle Physics and Astronomy Research Council. 



- 11- 

REFERENCES 

Antia H. M. & Chitre S. M., 1995, ApJ, 442, 434 

Bahcall J. N., 1988, Neutrino Astrophysics, Cambridge University Press, Cambridge. 

Barnett V. & Lewis T., 1984, in Outliers in Statistical Data, John Wiley & Sons, Chichester, 
p. 167 

Basu S., Christensen-Dalsgaard J., Schou J., Thompson M. J. & Tomczyk S., 1996, Bull. 
Astron. Soc. India, 24, 147 

Book S. A., 1978, in Essentials of Statistics, McGraw-Hill, London, p. 85 

Chaplin W. J., Elsworth Y., Isaak G. R., Lines R., McLeod C. R, Miller B. A. & New R., 
1996, MNRAS, 282, L15 

Christensen-Dalsgaard J., 1989, MNRAS, 239, 977 

Christensen-Dalsgaard J., 1996, in Proceedings of Workshop on Theoretical & 

Phenomenological Aspects of Underground Physics (TAUP '95), ed. Fatas M., Nucl. 
Phys. B. Suppl., in press 

Christensen-Dalsgaard J. & Dappen W., 1992, Astron. Ap. Rev., 4, 267 

Christensen-Dalsgaard J., Proffitt C. R. & Thompson M. J., 1993, ApJ, 403, L75 

Christensen-Dalsgaard J., Dappen W., et al., 1996, Science, 272, 1286 

Cox A. N., Guzik J. A. & Raby S., 1990, ApJ, 353, 698 

Dixon W. J., 1951, Ann. Math. Statist., 22, 68 

Dziembowski W. A. & Goode P. R., 1996, A&A, in press 

Elsworth Y., Howe R., Isaak G. R., McLeod C. P. & New R., 1990a, Nature, 347, 536 

Elsworth Y., Howe R., Isaak G. R., McLeod C. P., New R., 1990b, Nature, 345, 322 



- 12- 

Elsworth Y., Howe R., Isaak G. R., McLeod C. R, Miller B. A., New R., Speake C. C. & 
Wheeler S. J., 1994, ApJ, 434, 801 

Elsworth Y., Howe R., Isaak G. R., McLeod C. P., Miller B. A., Wheeler S. J. & New R., 
1995, in: GONG '94: Helio- and Asteroseismology From the Earth and Space, eds. 
Ulrich R. K., Rhodes Jr. E. J. & Dappen W., ASP Conference Series, vol 76, p. 76 

Evans D. J. & Roberts B., 1992, Nature, 355, 230 

Iglesias C. A., Rogers F. J. & Wilson B. G., 1992, ApJ, 397, 717 

Libbrecht K. G. & Woodard M. F., 1991, Science, 253, 152 

Regulo C., Jimenez A., Palle P. L., Hernandez F. P. & Cortes T. R., 1994, ApJ, 434, 384 

Rogers F. J., Swensen F. J. & Iglesias C. A., 1996, ApJ, 456, 902 

Tassoul M., 1980, ApJS, 43, 469 

Topping J., 1962, in Errors of observation and their treatment. The Institute of Physics & 
The Physical Society, Chapman & Hall Ltd., London, p. 88 



This manuscript was prepared with the AAS L^T[t;X macros v4.0. 



-13- 



Table 1. Comparison between BiSON and model fine-structure data 



Model Co (^Hz)"" ci (^Hz)'' Mean residual (a^Hz) 



do{n) 



Elsworth et al. 9.00 ± 0.06 -0.290 ±0.030 

Data, this paper 9.017 ±0.023 -0.3255 ± 0.0073 

Al 9.266 ±0.019 -0.3001 ± 0.0050 -0.170 ±0.013 

B^ 9.033 ±0.019 -0.2980 ±0.0053 0.050 ± 0.016 

C^ 9.143 ±0.019 -0.2953 ±0.0051 -0.081 ± 0.013 



rfi(n) 



Elsworth et al. 15.60 ±0.10 -0.460 ±0.060 

Data, this paper 15.903 ± 0.017 -0.4004 ± 0.0071 

A^ 16.103 ±0.020 -0.4165 ±0.0053 -0.201 ± 0.023 

B^ 15.789 ±0.015 -0.4112 ± 0.0041 0.120 ±0.024 

C^ 15.951 ±0.017 -0.4166 ±0.0045 -0.019 ± 0.024 

Tine-structure data fitted to equation di{n) = cq + ci • (n — 21), over range 
15 < n < 27 in order to facilitate comparison with Elsworth et al. (1990a). 

'^Weighted mean difference residual in sense BiSON minus model data. For 
do("-), nican calculated over range 9 < n < 28; for di{n), over range 11 < n < 27. 

^A "standard" solar model, with no settling effects. This is model 1 from 
Christensen-Dalsgaard, Proffitt & Thompson (1993). 

^Model 4 from Christensen-Dalsgaard, Proffitt & Thompson, incorporating 
settling with turbulent mixing. 

^Reference model 'S' from Christensen-Dalsgaard, Dappen et al. (1996), which 
incorporates settling effects, but neglects turbulence. 
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Fig. 1. — The upper [(io(^)] and lower [(ii(?2)] plots show residuals generated by subtracting 
the appropriate best straight-line fit cq + Ci ■ (n — 21) from the mean fine-structure spacings. 
The solid lines on each plot show the residuals from a "standard" solar model (A), while the 
dot-dashed (B) and dotted (C) lines show the residuals from models which include helium 
and heavy element settling (see § ^. 



